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ABSTRACT 

We study the photometric properties of stars in the data archive of the Sloan Digital 
Sky Survey (SDSS), the prime aim being to understand the photometric calibration over 
the entire data set. It is confirmed that the photometric calibration for point sources 
has been made overall tightly against the SDSS standard stars. We have also confirmed 
that photometric synthesis of the SDSS spectrophotometric data gives broad band fluxes 
that agree with broad band photometry with errors no more than 0.04 mag and little 
tilt along the wide range of colours, verifying that the response functions of the SDSS 
2.5 m telescope system are well characterised. We locate stars in the SDSS photometric 
system, so that stars can roughly be classified into spectral classes from the colour 
information. We show how metallicity and surface gravity affect colours, and that stars 
contained in the SDSS general catalogue, plotted in colour space, show the distribution 
that matches well with what is anticipated from the variations of metallicity and surface 
gravity. The colour-colour plots are perfectly consistent among the three samples, stars 
in the SDSS general catalogue, SDSS standard stars and spectrophotometric stars of 
Gunn & Stryker, especially when some considerations are taken into account of the 
differences (primarily metallicity) of the samples. We show that the g — r - inverse 
temperature relation is tight and can be used as a good estimator of the effective 
temperature of stars over a fairly wide range of effective temperatures. We also confirm 
that the colours of G2V stars in the SDSS photometric system match well with the Sun. 

Subject headings: 



1. Introduction 



The Sloan Digital Sky Survey (jYork et al.ll2000i . SDSS) provides a large data base of photome- 
try for galaxies and stars, over 350 million objects in the 11663 square degree survey area, that can 
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be used for a great variety of astrophysical science (jAbazaiian et al.l 120091 ) . Besides its scale, an- 
other very important feature is the special care that has been paid to maintain the accuracy of the 
data, especially their photometric accuracy, from instrumentation throughout the eight year period 
of the survey. The photometric system is different from systems which have been traditionally used, 
so it must be internally defined consistently over that period. In addition, the survey attempted 
to apply photometric calibration to the over 1 million spectra obtained during the survey, and we 
will study this calibration as well. 

One of the important issues is the consistency between spectrophotometry and broad-band 
photometry with the five passbands, in the sense that the integral of spectrophotometric data 
with the determined response functions should give the proper broad-band flux. This has been 
a perennial problem inherent with the data bases available to date. Synthetic colours obtained 
by integrating spectrophotometric data often give systematically increasing error with colours as 
obtained by broad band photometry, typi cally of the order of 0.1 mag, or sometim es larger, for stars 
for a range in B — V from to 1 (see e.g.. iFukugita. Shimasaku fc Ichikawalll995l ). This is probably 
largely to be attributed to a poor characterisation of the response functions of conventional detector 
systems. 

In order to minimise these systematic errors and make the broad-band photometry as well 
specified as poss ible, much effort was invested in the SDSS to characterise the response functions 



(|Doi et al.ll2010r ). in addition to a chain of stellar calibrations (I Tucker et al 



2006 



Stoughton et al 



2002). The response functions for SDSS photometry have been measured on site a number of times 
during the survey, which e nables us t o cha racterise their seasonal and secular variations. This 
work has been published in iDoi et al.l (|2010i ) , where the representa tive response functions are also 



presented as the Reference response functions of the 2.5m telescope (jGunn et al. 



1998 



20061 ). These 



average response func tions are r e comm ended to be used unless the user's interests are in fine details. 
It is demonstrated in lDoi et al.l (120101 ) that the stellar data taken in the early period of the survey 
are consistent with those acquired in the later period within 0.01 mag after calibration, including 
the data for the u passband which has suffered from a significant secular change in the system 
response and, therefore, the expected variations in photometry have been well compensated by the 
calibration process. It is concluded that the errors due to time varying changes, both seasonal and 
secular, in the response functions and calibration procedures that use several different systems are 
no more than 0.01 mag as a whole. 

The other special feature of the SDSS observation is hourly monitoring of the atmospheric 
extinction with the ancillary 50-cm Photometric Telescope (PT). A set of the SDSS standard stars 
were observed every night, together with hourly observation of stars in the transfer fields from PT 
to the 2.5-m main telescope. The SDSS standard stars are chosen as uniformly as possible both 
over the northern sky and in colour space s o that they are not inclined to give weight towards 



photometry, see lHogg et al 



particular class es of stars (Smith et al.l 120021). For oth e r efforts to mainta in the accuracy of the 



(|200lh : lTucker et al.l (l200fih : llvezic et al.l J^Km]) 
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In this paper we study the photometric accuracies and properties of the data presented in the 
SDSS catalogues, which are already in the public domain. The purpose of this paper is to examine, 
in particular, (i) the consistency between spectrophotometry and broad-band photometry, including 
the verification of the response functions for broad band photometry, (ii) the overall accuracy of 
photometric calibrations against standard stars, (iii) the location of stars in colour space in the SDSS 
photometric system and its variation with the atmospheric parameters of stars, such as effective 
temperature, surface gravity and metallicity, (iv) colours of galaxies relative to stars, and (v) 
brightness and colour of the Sun in the SDSS photometric system. One of the important purposes 
of (iii) is to locate stars of known spectral type in colour space so that approximate spectral type 
can be inferred from colours of the S DSS photometric s ystem. We also give effective temperature for 
those stars. The reader is referred to lLenz et al.l (|1998l ) for an earlier study along these lines. There 
is also a substantial amount of work done for observations of stars and for determinations of the 
atmospheric parameters in Sloan Extension fo r Galactic Understanding and Exploration (S EGUE) 
under the SDSS project fsee. lLee et all f|2008al Jbh : IXllende Prieto et all (|2008h : lYannv et all ([20091 )). 
Here, the problem is that the SDSS acquired spectra of stars fainter than g ~ 14, whereas the stars 
for which high resolution spectra are available needed to determine the atmospheric parameters 
are brighter, such as V < 14 (NB: V ~ g — 0.1): hence, the overlap of the samples is small, which 
makes well-calibrated stellar science with the SDSS not so straightforward. This also remains as 
our problem and we attempt to circumvent it by the use of bright SDSS standard stars, or of 
spectrophotometric standard stars after verifying their brightness in the SDSS photometric system. 



The procedure of SPSS broad band photometry is described in iTucker et al.l ( 200a ): see also 
Stoughton et al.1 (J2002J). In brief it is based on the network of 158 standard sta rs spanning a wide 
range of colours, which were measured by separate observations made at US NO (ISmith et al. 12002) 



using filters and a detector system pseudo-equivalent to the SDSS system (iDoi et al.ll20l"of ) with 



the colour tansformation between the two systems taken into account to offset the slight difference 
in their response functions. The zero point is ultimately set by the subdwarf BD+17°4708 whose 
brightness is defined by a synthetic calculati on of the spectrophoto metric data ()Oke h Gunnlll983l ; 
Okel 1 19901 ) in the AB95 magnitude system (jFukugita et al.lll996l ). The main telescope data are 
calibrated using the set of stars which are observed simultaneously with both the main telescope 
and an ancillary 50 cm Photometric Telescope, the latter being also used to observe the SDSS 
standard stars whose brightness is set by the work at USNO, and to monitor hourly atmospheric 
extinction during survey operation. 

Spectroscopy was carried out with two fibre-fed double spectrographs covering the wavelength 
range 3800 to 9200A with two 2048 2 -pixel CCD detectors in each spectrograph. The spectro- 



1 The full data rele ase (DR7) uses a global reca libration of the photometry using all the redundancies in the sample, 
called ubercalibrationl Padmanabhan et al-l lioOSl. We use here the standard pre-DR7 calibration; the mean difference 
between photometry with the two calibrations is of the order of 0.001 mag and the distributions are well described by 
Gaussian distributions with standard deviation 0.01 mag (g,r,i passband) and 0.02 mag (it and z passbands); These 
differences are smaller than we are concerned with in this paper. 
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graphs are fed with 640 fibers plugged into an alu minum plate in t he focal surface, positioned with 
the astrometry obtained from the imaging data ([Pier et al.ll2003l ). The mean resolving power is 
AA/A ~ 1850 — 2200 depending on the wavelength. Spectrophotometric calibration is carried out 
using spectrophotom etric standard stars sel ected in each field. Standard stars were intended to be 
typically F8 d warfs (Abazaj ian et alj 120041 bu t the set in actuality contains a variety of similar 
spectral types (lAdelman-McCarthy et al.ll2008l ): the spectra are fit to models given by Gray et al 
(2001), and the model spectra are used for calibration. Spectrophotometric calibration also requires 
match between spectrophotometric fluxes and the broad-band photometr ic PSF fluxes. This is done 



with a single calibration determined for each camera in a given plate (jAdelman-McCarthv et al 
20081 ). 



The data acquired in t h e SPSS projects are all published in 7 da t a releases, DR1— DR7 



(jAbazaiian et all 1200.1 [2004, 



2005 



Adelman-McCarthy et all 12009 . 120071 . 12008 ; I Abazajian et al 



20091 ). All DR data are cumulative, and DR7 stands for the final data release of the SDSS I and 
II, which employ the original instrumentation of the SDSS. We use the data set from DR6 in the 
study presented in this paper. The photometric calibration is made against the set of standard 
stars. Reddening corrections are not applied for t he stars. The stars are located mostly in low 
extinction regions. If we take the extinction map of lSchlegel et al.l ()1998l 1. the peak of the distribu- 
tion is at (B — V) ~ 0.02 with the mean 0.04 and the variance 0.03; only 2.6% of stars are located 
in the region E(B — V) > 0.1. Note that this reddening is the maximum value along the line of 
sight and the actual reddening for objects within the disk of the galaxy is smaller. 



2. Comparison between spectrophotometry and broad-band photometry 

In this paper we consider synthetic fluxes, obtained by integrating with the d\nv measure over 
the spectrophotometric data with the response functions for each filter. The flux calibration of 
the spectrophotometry was done by comparing the spectrophotometric data with the PSF fluxes 
of the standard stars in each spectroscopic field (actually two fields independently, one for each 
spectrograph.) So, it is important to confirm that synthetic spectrophotometric fluxes indeed agree 
with broad band photometric fluxes within tolerable errors for the majority of stars across the wide 
spectral range. 

In Figure [U we present the difference of synthetic colours and photometric colours for (a) g — r 
and (b) r — i as a function of colours of stars. Stars are taken from the observation for southern 
equatorial stripe (stripe 82), and the psf magnitude is adopted for broad band photometry. The 
error bar shows the variance in each colour bin, 0.2 mag in g — r or r — i. The offset of the mean 
is < 0.05 — 0.07 mag for most of the colour bins, and it increases to 0.08 mag for red stars with 
g — r ~ 1.3. The rms scatter we have seen is about 0.1 mag, and it goes up to > 0.15 mag towards 
redder stars with g — r > 0.7. This is compared with the rms error quoted for the spectrophotometric 
calibration, 0.04 mag (DR7). The smaller scatter for bluer stars offsets the larger rms scatter in 
longer wavelengths for the total sample. 
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The u and z passbands are out of the range of our spectroscopy, which spans from 3800A to 
9200A and therefore are not shown. We adopt the SDSS spect roscopic d a ta an d use the refer- 
ence response function for the SDSS 2.5 m telescope presented in iDoi et al.l (|2010l ) referred to the 
standard 1.3 airmass value for the atmospheric extinction corrections. 

The global tilt of the offset across colour is small. After 3<7 rejection of highly deviant data 
points, the mean gives A(g - r) spoc _p hot o = 0.026(5 - r) + 0.008 and A(r - i) sp ec- P hoto = 0.026(r - 
i) — 0.018 for the range —0.4 < g — r < 1.8 and —0.4 < r — i < 1.8, respectively. This implies that 
the response curves are well characterised. Synthetic colours may be substituted for broad band 
measurements at least in the mean over the stars in the sample and for non- variable stars. We 
consider that this is an improvement, crudely speaking, of a factor of 5 in accurac y compared with 
previous 'best' spectrophotometry-photometry data sets (e.g. jFukugita et al. Highly devian t 

data points are likely to be variable stars or double stars, as found for subsamples (jSesar et al.ll2007l ). 



3. Colours of stars in the SDSS data base 

In Figure [2] we plot u — g of 290647 stars with 16.95 < r < 17.05 against g — r taken from 
the full DR6 catalogue. Regions of densely populated points are represented by contours, which 
correspond to an increase of density by \/T0 in each step. Those denoted by individual points are 
2660 stars, or 0.9% of all stars. The catalogue contains an additional 583 stars that fall outside 
th e range and are no t plotted. We also superimpose the data o f 158 SDSS stan d ard s tars given 



m 



Smith et al.l (J2002J) and 175 spectrophotometric stars given bv lGunn k, Strvkerl (|1983l . hereafter 



GS83). The colour of GS83 stars is derived by synthetic calculations with the reference response 
function of the 2.5m telescope system at 1.3 airmass after taking the redleak of the u- filter into 
account^!. Since the GS stars are mostly located in the low latitude field and receive large reddening, 
dereddening correction is applied as in GS93. It is noted that the reddening model used in that 
paper is probably oversimple, although the resulting colour does not differ too significantly from 
what is calculated using a more standard extinction curve. 

There are 5 GS83 stars that fall beyond the figure frames and are dropped. All SDSS standard 
stars are shown. The redleak affects u — g colour by no more than 0.02 mag for stars with g — r < 1, 
while it increases to 0.2 mag close to the red edge of the contour, g — r ~ 1.5. Recognition of this 
correction improves agreement between the GS83 stars and other stars, if only slightly. The three 
solid squares around g — r ~ 0.28 represent the metal poor F subdwarfs, BD+17°4708, BD+26°2606 



Lenz et all (|1998| ) have also given SDSS colours for the GS83 stars using synthetic calculations. Their colours, 



however, do not agree well with ours, especially towards red stars (the difference becomes as much as 0.5 mag for 
u — g colo ur with q — r ~ 1 star s). This is much more than we expect fr om the difference between the two response 



functions, iFukugita et al. I |l996h , which was adopted by Lenz et al., and lDoi et all I |2010l ) adopted here. We expect 
that the two lead to a difference generally smaller than 0.02 mag and, at largest, < 0.04 mag for u — g colour of red 
stars. We could not find the origin of this gross discrepancy. 
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and BD+21°0607, the latter two of which were added to the firs t to set up the network of F- 
subdwarf standards approximately 8 ho urs apart rtSmith etalj 120021). Accurate spectrophotometry 
is available only for the first two stars (|Oke fc Gunnlll983l ; IOkd ll990). but the third is sufficiently 
similar to the first two that the adopted linear interpolation to obtain its fluxes should be adequate. 

We confirmed that the contours of the distribution of stars change only moderately when we 
take stars of r ~ 18 or 19 instead of r ~ 17. The distributions change in detail, doubtles s due to the 



changes in stellar popu lations from old disk to thick disk and halo as one goes fainter (jChen et al 



200ll ; llvezic et al.ll2008l ). The contours, however, overlap very well for stars bluer than g — r ~ 1 for 
the 3 samples that differ in brightness, indicating that the calibration is robust against brightness 
of stars. 

We see that both SDSS standard stars and GS93 stars are located near the centre of gravity 
of the distribution of stars in the SDSS general sample for g — r < 0.8. This also serves, albeit 
indirectly, as the evidence that spectrophotometry of the GS83 stars and spectroscopic synthesis 
give correct broad band colours to the level we have required. We see for red stars with g — r > 0.7 
some shift of the loci of these stars towards the upper side, u — g being redder by ~ 0.3 mag. On 
the other hand there appears a second locus of the GS83 stars, which is shifted downwards by ~ 0.1 
mag with respect to the centre of gravity. We ascribe this to widening of the gap in colour space 
between the distribution of luminosity class III and V stars for g — r > 1, as shown in the next 
section. For red stars we observe that the locus of the GS83 stars is more heavily populated around 
the upper side, to which stars are driven by high metallicity (nearly solar) of the GS83 stars. In 
general, the scatter becomes larger towards the redder end, as a natural wilder variation of spectra 
due to molecular bands developed in later K to M stars for g — r > 1. The span of the u — g colour 
range against the variation of metallicity, and also of gravity, increases towards the reddest stars. 

For very blue stars with g — r < —0.2, the locus of the GS83 stars does not agree with that 
of SDSS standard stars and that of the general SDSS stars, by 0.1 mag in g — r or as much as 
0.4 mag in u — g. This is significantly larger than could arise from the different response functions 
of SDSS specifications. One might suspect that this displacement is due to reddening, especially 
given the crudity of the GS reddening corrections. In fact, however, the displacement corresponds 
approximately to E(B — V) ~ 0.15. Such large extinction is very rare in the SDSS data, but is 
present in essentially all the blue GS stars, and is poorly determined. However, if such corrections 
are systematic and are made, we would expect shifts in i — z vs g — r, for which the shift (see Figure 
[6] belo w), if any, takes p lace in a way opposite to that expected. Note that the standard extinction 
curve (| O 'Donnelll 1 1 994! ) gives reddening for r — i and i — z smaller than that for u — g only by a 
factor of two (see arrows for extinction in Figure [8lfl0l below) . In r — i vs g — r the two loci overlap 
at a good precision; see Figure [5] below, but the extinction for this colour is nearly parallel to the 
stellar locus, and it is not easy to separate out extinction. We thus do not find an easy explanation 
that accounts for the shift in u — g vs g — r plane between the GS83 stars relative to the SDSS 
stars. We do note, however, that the populations of the blue GS stars and the blue SDSS stars are 
completely different; the blue GS stars are all main-sequence disk Population I objects, and the 
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blue SDSS stars are essentially all low-mass, hot, highly evolved Population II objects of a variety 
of types, so the discrepancy may not be surprising. We note also that the locus of the SDSS stars 
overlaps well that of the SDSS standard stars. 

To scrutinize the possible offset among the populations in three different star samples, i.e., 
the sample in the SDSS general catalogue, the SDSS standard star sample and the GS83 sample, 
we plot in Figure [3] the distribution of the three populations, projected onto the line orthogonal 
to (u — g) = 2.15(g — r) + 0.26 which follows the locus of the stellar colour for 0.2 < g — r < 0.8, 
indicated as the line segment in Figure El The distribution of the three F subdwarfs is also shown. 
The distance along this base line measured from the locus of the population is denoted by d ug (mag) 
with the positive direction towards the lower right (increasing g — r and decreasing u — g). We 
see that all populations are distributed very well with respect to the common zero point, although 
the distribution for the SDSS standard stars is somewhat skewed towards the bluer side, which can 
be understood by a larger weight of lower metallicity stars in that sample. No alarming offset is 
visible in the zero point among the four samples beyond ~ 0.05 mag in d ug . 

The distance d ug is correlated with metallicity. 21 out of the 158 SDSS standard stars and 
44 out of the 175 GS83 stars are given esti mates for metallicity with high resolution spectroscopic 



observations (jCayrel de Strobel et al.ll200ll ). Figure [4] shows those stars having colours in the range 
0.2 < g — r < 0.8 plotted in the plane of metallicity versus d ug , indicating the correlation of 
metallicity with d ug . It shows that the GS stars plotted in this figure have metallicity close to 
solar, while many of the SDSS standard stars are indeed metal poor [Fe/H|< —1, as expected from 



the fact that they are mostly high latitude halo objects (llvezic et al.ll2008l ). Although these are 
partial samples of the GS83 and the SDSS standard stars, this explains the difference of the two 
samples we saw in the u — g vs. g — r plane. We also plot the bundle of lin e segments which show 



the prediction of the Kurucz atmosphere model (jCastelli Kuruczl |2004| ) for stars with F2V to 



G8V, metallicity from [Fe/H]=0.5 being at the top to —1.5 at the bottom. 

A similar plot is shown in Figure [5] for r — i vs. g — r. The distribution is narrower and the 
agreement among the three distributions is much tighter than we have seen for u — g vs. g — r. The 
SDSS standard stars are located very close to the centre of gravity of the general star sample. The 
GS83 stars are also located close to the centre of gravity: we observe a slight displacement between 
the SDSS standard star sample and the GS83 star sample, which is of the order of the difference 
we expect from the different metallicity. The three F subdwarfs are also on the top of the loci. The 
difference in the loci between the luminosity class III and V is negligibly small, as we see in the 
next section more quantitatively. 

A similar observation holds for the plot of i — z against g — r shown in Figure El We do not 
observe displacements among the three populations, again except for some disagreement for very 
blue stars. Otherwise, the distributions of the standard stars and GS83 stars agree well with those 
of the stars in the SDSS general catalogue. 



These considerations mean that the stellar loci are well defined in the SDSS photometry and 
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depend little on the sample taken, and also verify that the photometric calibration has been made 
well in colour against the standard stars. This does not mean that the broad band flux is accurate 
in physical units (AB magnitude system), however. 



4. Variations of colours and stellar atmosphere 



The main parameters that regulate colours of stars are effective temperature T e g, surface 
gravity g and met allicity [Fe/H]. We first show that T e g is well represented by g — r colour (see also 
Ivezic et al.l (|2008l )). The most accurate measurement of T e g- is made with Michelson interferometry 
of stars to measure their diameters. This has been done, however, only for very bright stars and the 
sample is very limited. The sample may be exten ded by stars for which tem peratures have been 
obtained from the infrared flux method (IRFM) (jBlackwell Shallia 119771 ). which is calibrated 
against stars with Michelson interferometry measurements of high accuracy. The IRFM suffers 
from less model dependence; temperatures from the IRFM method are supposed to be accurate 
at the 1% level. The stars with effective temperatures estimated by IRFM are still too bright 
for the main sample of SDSS, but the GS83 sample has stars whose effective temperature has 
been estimated with IRFM. The argument in the previous sections ensures that photometry from 
spectrophotometric synthesis for GS83 stars is sufficiently accurate. 

Figure [7] shows col ours of 15 GS 8 3 sta rs against inverse effective temperature that is given in 
the IRFM catalogue of Idi Benedetto! (jl998l ). Among the panels that show u — g, g — r, r — i and 
i — z colours the correlation with g — r colour is tightest: the stars are well fitted with the simple 
formula 

i hi i 

(1) 



T eff /10 4 K 

with an rms scatter around the fit a, 



1.09 



g-r + 1.47 

g -r — 0.04 or cr(T e ff) = 93K. Note that the star sample 
contains different luminosity classes, III (7 stars), IV (4 stars) and V (4 stars), and metallicity 
is not differentiated. We do not observe any systematic difference for g — r among the different 
luminosity classes. We note, however, that this colour-temperature relation is empirically verified 
only in the range, 1.25 < 10 4 K/T eff < 2.6 or 3850K < T cS < 8000K. The use of the Castelli-Kurucz 
atmosphere model (indicated by the dotted curve in the second panel of Figure [7|) suggests that 
this inverse linear relation may break down for T e g- > 9500K: for those stars with temperature 
higher than T e g ~ 8500K, T e g is higher than eq.([l]) indicates for given g — r. A departure from 
eq.(H]) is also indicated towards lower temperature below 3850K. 

The expre s sion w e derived above agrees well with the log T e g- — (g — r) colour relation given 
by llvezic et al.l (|2008h but only in a finite segment, -0.3 < g-r < 0.7, or 5000 < T cff < 8000K. 
The latter relatio n gives temperatu re significantly lower than the IRFM estimate for stars redder 
than g — r ~ 0.7. iLee et al.1 (j2008al ) use a polynomial which is third order in g — t to estimate Teff 
from g — r in the SEGUE stellar parameter pipeline. It agrees with eq. (PQ) within 0.02 dex for the 
range -0.6 < g - r < 1.2 (4200K < T cff < 12500K). 
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We see some systematic shifts among different luminosity classes for the u — g vs. inverse 
temperature plot; e.g., all 4 class V stars lie below the fit line, and 6 of the 7 class III stars lie 
above the line. The i — z colour is also well correlated with temperature, the dispersion of i — z 
being also small (0.03). The corresponding dispersion of temperature, <r(T e ff) = 190K, however, is 
twice as large as that for the g — r inverse temperature relation due to the shallower slope. 



Metallicity has been estimated with high resolution spectroscopy (jCavrel de Strobel et al 



200 ll ) for 9 among the 15 stars given in this figure and it ranges from [Fe/H]=— 0.69 to 0.31. 
We do not observe a systematic trend with metallicity in the g — r vs. T e g plot. The variation 
with metallicity is significantly smaller than is seen with B — V, since the metal lines are fewer 
and weaker in the g (and r) passband than in the B passband. The scatter around the colour 
temperature relation is nearly as small as the one that would be obtained after the metallicity 
effect taken into account in the B — V - T P ff relation, or the V — K - T P ff relation for which the 



smallest scatter has been claimed ([Cohen et al.lll978l ; Idi Benedettolll998l ). Some metallicity trends 



are observed for the u — g - T e g plot in the way one would expect, but the present sample is too 
scanty to derive the metallicity dependence in a meaningful way. We conclude that g — r colour is 
a good indicator for (inverse) temperature, much better than B — V because of weak metal lines 
in the g and r passbands. 

We repeat in Figure El the u — g vs. g — r colour-colour distribution shown in Figure EJ In 
panel (a) the loci of GS83 stars are represented by the two curves, the solid curve for luminosity 
class V and the dashed curve for class III. These curves are drawn by spline-interpolating the 
mean positions of the GS83 stars that are classified as III and V. The curve segments that are not 
constrained well by the data are skipped in the figure. The ellipses are the areas where GS83 stars 
with some specified types, AO and Al stars (denoted as AO— 1), F0 and Fl stars, GO, Gl and G2 
stars, K0 stars and M0 stars fall (as indicated), where III and V are not distinguished. The size of 
the ellipses corresponds to the sample variance. The two curves cross at around g — r = 0, but the 
data are scanty for class III and this crossover is yet to be examined with a larger data set. We see 
that the gap between the two curves widens for g — r > 1, which is likely to accoun t for the gap 



betwe en the two branches in the GS83 sample as we saw earlier, and as also noted by lYanny et al 



(|2009l ). We indicate at the bottom right of the panels how reddening affects co lour. The length 



of the arrow corresponds to E{B — V) = 0.3, and the extinction curve we used (| Q ' Donnelll 1 1 9941 ) 
gives k\ = A X /E(B - V): k u = 5.155, k g = 3.793 and k r = 2.751. 

We give in Table Q] the effective temperature and colours for given spectral types for stars in 
the GS83 sample, where the use is made of our eq. (1) to estimate the effec tive temperature o r, 



where available, the temperature is adopted from the IRFM measurements (Idi Benedetto! 11998). 
We calculate the mean over the available stars, although for most listings only one or two stars are 
available in the GS83 sample. 

It is expected that colour is affected by metallicity. To show the range in u — g vs. g — r 
colour space that is covered by the variation of metallicity, we indicate it in Figure E^b) with 
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line segments where the two edges correspond to [Fe/H] =+0-5 (upper edge ) and — 1 (lower edge) 
for luminosity class V, taking the atmosphere model of ICastelli fc Kuruczl feooi fl The asterisk 



on the segment shows the position of [Fe/H]=+0. These line segments are shifted for luminosity 
class III stars, somewhat irregularly for K and M stars. We expect that the Kurucz atmosphere 
models may give the relative metallicity-dependent position with reasonable accuracy, though the 
absolute position may somewhat be shifted. The models properly represent the range covered by 
the variation of metallicity. This shows that the width of the distribution of stars in the general 
catalogue is consistent with the variation covered by the range of metallicity —1 <[Fe/H]< along 
with further widening of the range due to the mixture of luminosity class III and V stars. 

A similar figure is given for the r — i vs. g — r colour distribution in Figured The splitting 
of tracks between luminosity classes III and V is very small: the two curves are nearly degenerate 
except for red stars with g — r > 1.2. The variation due to metallicity is also reduced although 
it is still substantial in g — r. This generally smaller variation makes the discrepancy between the 
observation and the Kurucz atmosphere model, indicated by line segments in Figure [9b, somewhat 
more apparent. The Kurucz atmosphere model fits the curve representing the data only at its 
lowest metallicity edge, although the discrepancy is only about 0.03 mag in r — i. The trend is 
similar for the i — z vs. g — r colour distribution (Figure [T0"|h We see a small gap appearing between 
the two luminosity classes which is somewhat more conspicuous than is seen in r — i vs. g — r. 
Metallicity induced colour variation is larger for i — z, even larger than for u — g, for red stars. 

All of the above indicates that the SDSS system is relatively close to a true AB system; the 



suspected corrections to an AB system of a few hundredths of a magnitude (lAbazajian et alJ l2009) 
are small enough not to have been seen in this analysis. 



5. Colours of other objects 



To give the idea as to the feature of SDSS colours, we show in Figure [TT] the distribution of 
white dwarfs, the sam ple of which is taken fro m the 9300 spectroscopically confirmed white dwarfs 
(16.75 < r < 17.25) of lEisenstein et al.l (120061 ). They are white dwarfs with temperature typically 
> 7000K. Cooler white dwarfs are degenerate with main sequence stars, and more with subdwarfs. 



Selec t ions other than co lours are needed to find candidates for cooler white dwarfs (jKilic et al 
20061 : lHarris et al.1 120061) . The l ocus o f quasars is also added (solid curve) taking the composite 



spectrum of IVanden Berk et al.1 (|200ll ) as the fiducial and by redshifting it to z < 2.5. The lowest 
redshifts shown depend on colour: they are for u — g, 0.05 for g — r, and 0.35 for r — i. The locus 
of Type la supernovae at the epoch of their B-band brightn ess maximum is also shown (dashed 
curve) using the model spectrum of SALT2 ([Guv et al.l 120071 ) that is moved to non-zero redshift, 



Castelli fc Kuruc3 1 2004h gives the prediction significantly different from Kurucj 1 19931 ) for red stars with 



g-r > 1 



in this figure due to the inclusion of H2O opacities and the revision of TiO lines. 
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z < 0.5. Some parts of the locii lie outside the frame of the figure. 



The figure shows the distribution of white dwarfs well separated from normal stars due to 
their UV excess in the u — g vs. g — r plane, in so far as hot w h ite dw arfs (approximately with 
T > 7000K) are concerned, as efficiently used in lEisenstein et al.1 (120061 ). and to some extent also 
in the i — z vs. g — r plane. A strong degeneracy with normal stars is seen in the r — i vs. g — r 
plane. White dwarfs are degenerate with the blue tip of the distribution of normal stars with these 
colours. 

Low redshift quasars, having u — g colors much bluer than most stars, are also well separated 
from normal stars on the u — g vs. g — r plane till quasars reach z ~ 2.5, where they mingle 
into the stellar locus. The confusion of quasars with white dwarfs may take place in the u — g vs. 
g — r plane, though quasars are somewhat redder in g — r: this degeneracy is lifted if the i — z 
vs. g — r colour is used in addition. This demonstrates the advantage of multicolour space in 
the photometric selection of low redshift quas ar candidates. Th e ph otometric target selec t ion fo r 
quasars is discussed at length in the papers of iFan et al.l (|200ll ) and iRichards et al.l (|2002l . l2004j ) , 
and we do not discuss this problem further. 

The figure shows that Type la supernovae can be identified, well separated from the stellar 
locus, given the homogeneity of colours of Type la supernovae. The fiducial colour of the SALT II 
template at zero redshift, u — g = 0.476, g — r = —0.170 and r — i = —0.673, is compared with that 
derived from the SDSS supernova sample, which reads u — g = 0.591±0. 237, q — r = —0.093 ±0.15 7 
and r—i = — 0.636±0.088 where the error represents the sample variance (jYasuda &: Fukugitall2010i ). 



The diagram shows that Type la supernova locus moves monotonically with redshift, indicating that 
the photometric redshift can be estimated with a good confidence, given good colors at maximum 
light. 

The three panels of Figure [12] s how colou r s of m orphologically classified galaxies at z ~ 0. 
Galaxy colours are taken from iFukugita et al.l (|2007l ) based on the visual classification of 2250 
galaxies with r < 16 in the northern equatorial stripe. Magnitudes are Petrosian. The error 
bars indicate the variance of the samples for E, SO, Sa, Sb, Sc and Im. Galaxies occupy a colour 
space locus narrower than stars, with colours roughly corresponding to F5— K2 stars in the range 
g - r = 0.2 - 0.8. 



6. Brightness and colour of the Sun 



The brightness of the Sun is taken as a basic unit in astrophysical work. The measurement of 
the brightne ss of the Sun , however, is notoriously difficult due to its extreme brightness and angular 
extent (e.g.. lHayeslll985l ). The measured broad band colour of the Sun has been variously reported 
from (B — V )(7) = 0.62 to 0.68 by various investigators. The choice often taken is (B — V) & = 0.65 
(lAllenlll973h . 
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We adopt the absolute spectrophotometry of the Sun measure d by SOLSPEC in th e Atmo- 
spheric Laboratory for Application and Science (ATLAS) 3 mission (jThuillier et al.l 12003 ) . the ab- 
solute flux of which has been calibrated against a black body standard. We present in Table [2] the 



synthetic broad band magnitudes of the Sun f or V, B — V, q 



u - 



a — r , r — i and i — z. We take 



the B and V response functions evaluated bv lAzusienis Straizvsl (|1969l ). which seem to give the 
minimal offset of synthetic brig htness against Johnson broad band ph otometry among a number 
of response functions we tested (jFukugita. Shimasaku k, Ichikawalll995l ). We also pres ent synthetic 



values of brightness of the Sun using another spectrophotometric table compiled by IColina et al 
( 19961 ) in order to indicate the size of possible systematic errors. 



The agreement among the different spec tral energy dist r ibutio n (SED) is within about 0.02 
mag for the V band brightness. The SED of iThuillier et al.l (|2003l ) integrated with the reference 
response function of the SDSS 2.5m telescope with 1.3 airmasses, gives V = —26. 73 mag which 
agrees well with —26.75 ± 0.06 mag of the summary value of photometry, given by lHavesI (|1985l ) 
and is also close to the value, —26.74, adopted in Allen's Astrophysical Quantities (1973). The use 
of the SED from earlier missions of ATLAS also gives V band brightness w ithin 0.02 mag of the 
ATLAS3 value. This is, in fact, the order of magnitude IThuillier et al.l (|2003l ) claimed as the error 
of calibration. 

This result leads u s to expect that solar SEP can also be used to calculate broad band colours. 
The B — V colour from IThuillier et al.1 (J2003J) is 0.62 ±0.01. While this looks somewhat bluer than 
the conventional value, we should remark that B — V colour estimated from the colour temperature 
rela tions for main sequence star s with the solar abundance also lie bluer, in the range 0.616—0.635 
(see lSekiguchi fe Fukugitall2000l ). The colour often quoted, B — V = 0.650, is rather close to colour 
of G4 stars than G2, and is suspected to be too red. 

The modern SED of the Sun gives a reasonable agreement among the different data sets for 
both absolute flux and colour. We calculate with a reasonable confidence g — r, u — g, r — i and i — z 
colour of the Sun as given in Table [2] in the AB95 magnitude system. The u — g and g — r colours 
thus obtained (see the cross symbols in Figure [8]) fall close to the later-type edge of the 1 a ellipse 
representing GO— G2 stars of the GS83 sample in Figure El [u — g) & = 1.21 and (g — r) & = 0.45 
are close to the mean of two G2V stars in the GS83 sample, 1.22 and 0.43, respectively. The same 
is true also for the redder passbands: (r — i) Q = 0.13 and (i — z) & = 0.02 are compared with 0.12 
and 0.03 of the GS83 stars. The Sun thus matches very well with G2 stars. The brightness in the g 
passband is —26.46. Our V and g magnitudes translate to absolute brightness My(Q) = 4.84 and 
Mo(©) = 5.11 taking the distance modulus —31.57. For other colours, see Table 2. 



7. Summary 



We have studied photometric properties of stars given in the data archive of the SDSS, the 
primary question being whether the photometric calibration was, overall, properly done. We found 
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that, over the entire lifetime of the survey, the photometric calibration for point sources has been 
made tightly against the SDSS standard stars and colours of the stars are well defined. We have 
not identified a sample of stars which are significantly and systematically deviant from the SDSS 
standard stars. We have also found that the synthesised colour from the Gunn Stryker spectropho- 
tometric sample represents very well the colour of stars in the SDSS general catalogue and vice 
versa. Photometric properties are perfectly consistent mutually with each other among the three 
samples we studied. It is also gratifying that the the synthesis of the SDSS spectrophotometric 
data with the use of the reference response function of the 2.5m telescope of the SDSS give broad 
band fluxes that agree with broad band photometry of SDSS with little tilt along colours. 

We have also given the fiducial colours and temperatures for stars empirically given spectral 
types, and show how metallicity and surface gravity affect colours. This enables us to infer spectral 
type of stars when SDSS colours is given. We show that g — r colour can be used as a good 
estimator of the effective temperature of stars. The distribution of stars in colour space matches 
well with what we expect from the variations of metallicity and surface gravity. We also present 
the brightness and colour of the Sun with synthetic calculations, which shows that colour of the 
Sun matches well with G2V stars in the SDSS photometric system. 

A problem we encountered in the present study is that the sample of stars with well known 
temperature and metallicity is still too scanty, and bright stars that have SDSS photometry are small 
in number. This is especially depressing when we require these two quantities at the same time. 
There is a pressing need to acquire photometry with the SDSS passbands for bright star samples 
which have good measured atmospheric parameters, (composition, temperature, and gravity) to 
connect more tightly SDSS photometry with physical stellar parameters. 
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Table 1. Temperature and colour of stars. 



Type 


colour 


T e & 
IRFM 


u-g 


g-r 


r — i 


Type 


colour 


IRFM 


u-g 


g-r 


r — i 


rJ2V 


(11195) 




—0.08 


—0.50 


—0.37 


B2111 


(11180) 




—0.10 


—0.50 


—0.35 


B7V 


(9631) 




0.49 


-0.34 


-0.25 


B7III 


(9832) 




0.49 


-0.36 


-0.26 


AOV 


8715 




0.98 


-0.22 


-0.19 


A0III 












A5V 


7656 




1.08 


-0.04 


-0.09 


A5III 


7812 


7997 


1.10 


-0.07 


-0,13 


FOV 












F0III 












F4V 


6831 




1.11 


0.13 


-0.11 


F4III 












F9V 


5912 


5902 


0.90 


0.37 


0.10 


F9III 












GOV 


6008 




1.04 


0.34 


0.09 


G0III 












G2V 


5747 




1.22 


0.43 


0.12 


G2III 












G5V 


5596 




1.29 


0.48 


0.14 


G5III 












G8V 


5217 


5155 


1.51 


0.62 


0.21 


G8III 


4969 


5026 


1.86 


0.72 


0.23 


KOV 


5610 




1.19 


0.47 


0.15 


K0III 


4731 




2.14 


0.84 


0.28 


K4V 


4600 




2.06 


0.90 


0.29 


K4III 


4017 


4050 


3.09 


1.24 


0.43 


K7V 


3879 


4011 


2.47 


1.34 


0.54 


K7III 












MOV 


(3793) 




2.44 


1.41 


0.66 


M0III 


3828 


3846 


3.16 


1.38 


0.72 


M2V 


(3700) 




2.47 


1.48 


0.92 


M2III 


(3585) 




3.10 


1.57 


1.28 


M5V 


(3415) 




2.29 


1.72 


2.12 


M5III 


(3576) 




2.60 


1.58 


1.89 



Note. - - T e fj (colour) is estimated from eq.(l) and the applicability is suspected for values given in 
parentheses; see the text. 
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Table 2. Brightness and colour of the Sun. 



SED Data 


V 


B - V 


9 


u- g 


g-r 


r — i 


i — z 


Thuillier et al. ( 2003) 


-26.734 


0.622 


-26.459 


1.213 


0.451 


0.127 


0.015 


Colina et al. (1996) 


-26.748 


0.644 


-26.453 


1.234 


0.479 


0.106 


0.020 




My 




M u 


M g 


M r 


Mi 


M z 




4.84 




6.32 


5.11 


4.66 


4.53 


4.52 



Note. — Synthetic broad band absolute brightness obtained using Thuillier et al. 
(2003) 's SED. V magnitude is in the Johnson system, giving V = 0.03 for a Lyr. 
ugriz magnitudes are in the AB95 system. 
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Fig. 1. — Difference in the brightness of broad band colour between the synthesis of the spectropho- 
tometric data and the direct broad band photometry, m(spectro) — m(photo), is plotted against 
colours of stars acquired in the southern equatorial stripe (Stripe 82) for g — r and r — i. The lines 
are linear fits with the coefficients given in the figure. 
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Fig. 2. — Three hundred thousand stars in u — g vs g — r colour space taken from SDSS general 
star catalogue. Densely populated regions are shown by contours, and one step corresponds to 
VTU increase in the density with the outermost contour for 75 stars/(0.1 mag) 2 . Open circles show 
SDSS standard stars, and solid circles are the GS83 stars. Those denoted by solid squares are three 
BD F-subdwarfs, taken to give the observational flux standard. The oblique line segment along the 
contour (from g — r = 0.2 to 0.8) indicates the approximate positions of the centre of gravity for 
the star distribution used to define d ug used in later figures. 
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Fig. 3. — Distribution of stars with 0.2 < g — r < 0.8 in different samples projected to the axis 
that is perpendicular to the distribution of densely populated region, which is around u — g = 
2.15(5 ~~ r ) + 0-26 in g — r vs u — g colour space, shown in Figure 2. The abscissa is d ug (mag). 
Open histogram shows the stars in the general SDSS catalogue, the histogram hatched from up-left 
to down-right (red) shows stars in the GS83 sample, and that hat c hed f rom up-right to down- left 
(green) shows the SDSS standard stars presented in ISmith et al.1 (120021 ) . The thick horizontally 
hatched histogram at around the zero point in the abscissa shows 3 F subdwarfs taken as the 
fundamental standard at the observation. The ordinate is in arbitrary units. 
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Fig. 4. — Correlation of [Fe/H] estimated in ICayrel de Strobel et al.1 (120011 ) with d ug , i.e., the 



distance from the line (u — g) = 2.15(g — r) + 0.26, which represents the centre of gravity of the most 
densely populated region between g — r = 0.2 and 0.8, as defined in the text. The sample is taken 
from the GS83 stars (solid circles) and the SDSS standard stars (open circl es) . The line segments are 
the metallicity dependence predicted by the Kurucz atmosphere model (ICastelli h Kuruca 12004) 
from F2 to G9 stars in the main sequence. 
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Fig. 5. — Same as Figure 2, but for r — i vs. g — r. 
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Fig. 6. — Same as Figure 2, but for i — z vs. g — r. 
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Fig. 7. — u — g, g — r, r — i and i — z colours plotted against inverse temperature in units of 10000K. 
The GS83 stars for which temperature is estimated using the IRFM are shown. The solid circles 
(red) show stars with the luminosity class V, open circles (blue) stars with luminosity class III and 
solid triangles (green) stars with luminosity class IV. Lines are linear fits to the data and their 
coefficients given in the figure where parentheses are dispersion of the fit in colour units. In the fit 
x = 10 4 K/T e ff and y is the corresponding colour. The dotted curve in the panel for g — r is the 
prediction of the Castelli-Kurucz atmosphere model for main sequence stars. 
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Fig. 8. — (a) Distribution of stars in the SDSS general catalogue in u — g vs. g — r plane (the 
same as Figure 2 but with isolated points representing individual stars suppressed). Two curves 
are empirical loci of main sequence stars obtained by interpolating the GS83 stars with known 
spectroscopic types separately for luminosity class V (solid curve, red) and III (dashed curve, 
blue). In panel (a) ellipses show the position of the GS83 stars with types given as indicated in 
the figure with the size representing the variance of the subsample. The plus symbol is colour of 
the Sun given in the text. In panel (b) [ see next page ] line segments indicate variations due 
to metallicity from [Fe/H]=— 1 to +1 with 0.5 unit steps (blobs), from downstairs up, using the 
Kurucz atmosphere model. Asterisks show the position for [Fe/H]=0. 
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Fig. 8. — (b) continued 
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Fig. 9. — (a) Same as Figure [8] but in r — i vs. g — r plane. 




Fig. 9. — (b) continued 
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Fig. 10. — (a) Same as Figure [8] but in i — z vs. g — r plane. 
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Fig. 10. — (b) continued 
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Fig. 11. — Distribution of stars in (a) u — g vs. g — r plane, (b) r — i vs. g — r plane, and (c) 
i — z vs. g — r plane (the same as Figure [2J [5] and respectively, with isolated points representing 
individual stars suppressed). Small dots show hot white dwarfs. The solid curve shows the locus of 
quasars having the composite spectrum with redshift ~ to 2.5, with blobs at the 0.5 interval in z: 
the lowest redshift is shown with the asterisk, which is z = for (a), z = 0.05 for (b) and z = 0.35 
for (c). The dashed curve is the locus of Type la supernovae redshifted to z < 0.5 towards red. 
The blobs on the supernova curve are at a 0.1 interval in redshift, but those at z ~ is out of the 
range of the figure and omitted in (b) and (c). The blobs at z = 0.1 and 0.2 are also omitted in 
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Fig. 12. — Distribution of stars in (a) u — g vs. g — r plane, (b) r — i vs. g — r plane, and (c) 
i — z vs. g — r plane, where colours of galaxies are overlaid. The error bars represent the variance 
of colours in the morphologically classified sample, from the right to the left E, SO, Sa, Sb, Sc, Sd 
and Im. 



